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The Korteweg-de Vries (KdV) equation is a non-linear wave equation that has played a funda- 
mental role in diverse branches of mathematical and theoretical physics. In the present paper, we 
consider its significance to cosmology. It is found that the KdV equation arises in a number of 
important scenarios, including inflationary cosmology, the cyclic universe, loop quantum cosmology 
and braneworld models. Analogies can be drawn between cosmic dynamics and the propagation of 
the solitonic wave solution to the equation, whereby quantities such as the speed and amplitude 
profile of the wave can be identified with cosmological parameters such as the spectral index of the 
density perturbation spectrum and the energy density of the universe. The unique mathematical 
properties of the Schwarzian derivative operator are important to the analysis. A connection with 
dark solitons in Bose-Einstein condensates is briefly discussed. 
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I. INTRODUCTION 



The Korteweg-de Vries (KdV) equation [l| is the com- 
pletely integrable, third-order, non-linear partial differ- 
ential equation (PDE): 



d t u + d 3 u H ud x u = 0, 

u 



(1) 



where u = u(x,t), d t = d/dt, d 3 = d 3 /dx 3 , etc., uq is a 
constant and (x, t) represent space and time coordinates, 
respectively. This equation was originally derived within 
the context of small-amplitude, non-linear water wave 
theory and it is well known that it admits a solitonic 
wave solution of the form 



u = u A 2 sech 2 [A(cc - X 2 t)/2], 



(2) 



where the constant A/2 represents the wavenumber of 
the soliton The KdV soliton is characterized by the 
property that its speed and amplitude are proportional 
to the square of the wavenumber. 

The KdV equation has played a central role in diverse 
branches of physics, including nonlinear optics, atomic 
and nuclear physics, Bose-Einstein condensates and as- 
trophysical plasmas (see Ref. @ for a review.) As far 
as we are aware, however, it has not been discussed pre- 
viously within a cosmological context (although see Ref. 
Q). This is perhaps not surprising, given that the KdV 
equation is a third-order PDE in two independent vari- 
ables, whereas the field equations for spatially isotropic 
universes are second-order, ordinary differential equa- 
tions (ODEs). 

One the other hand, Eq. (JlJ can be reduced to the 
non-linear ODE: 



■ AV + u" H uu 

u 



— uu' = 0, 



(3) 

where a prime denotes d/d<f> and (f> = x — X 2 t repre- 
sents a 'wave-like' independent variable. The purpose 



of the present paper is to show that wave solutions to 
Eq. and in particular the soliton solution ([2]), arise 
in a number of cosmological settings, including the in- 
flationary paradigm, the cyclic universe scenario, loop 
quantum cosmology and braneworld models. We con- 
sider scenarios where the universe is dominated by a sin- 
gle, self-interacting scalar field, <f>. Interpreting the value 
of the scalar field in terms of a wave-like coordinate of 
the KdV equation then allows for a direct analogy to be 
drawn between cosmological dynamics on the one hand 
and solitonic behaviour on the other. Cosmological pa- 
rameters can then be identified with quantities such as 
the speed and amplitude of the corresponding wave. 

The paper is organised as follows. In Section 2, we 
discuss a connection between the KdV equation and the 
Schwarzian derivative operator that we employ in later 
Sections. We proceed in Section 3 to consider the classes 
of inflationary and cyclic universes that generate density 
perturbation spectra with a constant spectral index. In 
Section 4, we find that the scaling solutions of various 
braneworld and loop quantum cosmological scenarios are 
analogous to the KdV soliton. We conclude with a dis- 
cussion in Section 5 on a connection with dark solitons 
in Bose-Einstein condensates. Unless otherwise stated, 
units are chosen such that H = c = 1 and the Planck 
mass is normalised to mp — y8n. 



II. THE KDV EQUATION AND THE 
SCHWARZIAN DERIVATIVE OPERATOR 

The KdV equation ([3]) admits an auto-Backlund trans- 
formation, whereby a solution u = u(d>) can be derived 
from a given seed solution u = u(4>) [4j. For the special 
case where the seed is the trivial solution u — 0, such a 
transformation reduces to the condition that a solution 
to Eq. © is given by 



u {\ 2 - Ay 2 ), 



(4) 



2 



where A 2 is a constant and the function y = y(<j)) is a 
particular solution to the first-order Riccati equation: 



For example, the solution 
A 



A 2 
4 



y=- tanh(A0/2) 



(5) 



(6) 



generates the KdV soliton $2$. 

Solutions to Eq. ([5]) satisfy a third-order ODE given 

by 



S[y(<t>)] 



3 fv" 



y< 2 V y' 



A 2 
" 2 ' 



(7) 



as may be verified by direct differentiation of Eq. ([5]). 
The left-hand side of Eq. is the Schwarzian deriva- 
tive operator (often referred to simply as the Schwarzian) 
Q . This operator exhibits a number of remarkable prop- 
erties, one of which we exploit in the present work: it 
is the unique combination of derivatives that is invari- 
ant under a homographic transformation corresponding 
to the group of fractional linear transformations. This 
follows since, for a function y = y (</>), the composition 
N(y) = [ln(y')]' = y"/y' transforms under the inver- 
sion y — > 1/y such that N(l/y) — N(y) — 2y'/y. Some 
straightforward algebra then implies that the operator 
S(y) = N' — N 2 /2 is invariant under inversion. More- 
over, since S(y) = S{my + n) for any m, n G 5i, S(y) is 
invariant under the full group of fractional linear trans- 
formations. 

This implies that if y (<f>) is a particular 'seed' solution 
to the differential equation S[y] = f(4>) for some function 
f(<p), the general solution to such an equation is given by 



2/(0) 



ay ■ 



cy ■ 



d' 



a b 
c d 



eSL(2,R), (8) 



where a, b, c, d are constants such that ad — be = 1. 

To summarize, therefore, if a particular solution to the 
Schwarzian equation ([7]) for constant A can be found, the 
general solution ([5]) can be written down immediately. 
Restricting the general solution to satisfy the first-order 
constraint (J5J) then generates a solution (Q| to the KdV 
equation ([3]). We employ this result in the following Sec- 
tions for a number of cosmological models. 



III. COSMOLOGICAL DENSITY 
PERTURBATIONS AND THE KDV SOLITON 

A. The Inflationary Universe 

The inflationary scenario remains the cornerstone of 
modern, early universe cosmology. Whilst there is cur- 
rently considerable interest in multiple-field versions of 



the paradigm, our aim in this Section is to revisit the 
simplest version of the scenario, namely inflation driven 
by a single, minimally coupled, slowly-rolling scalar in- 
fiaton field, (f>. (For reviews, see, e.g., Ref. jaulO We 
consider the general class of models that generate a scale- 
invariant density perturbation spectral index (hereafter 
referred to as the spectral index) to lowest-order in the 
slow-roll approximation. Our main aim is to highlight 
some interesting mathematical features of the underly- 
ing differential equations that have not been previously 
discussed. 

The cosmological Friedmann equations in Hamilton- 
Jacobi form are given by 



3iT - 2H' 



v (<!>), ^ 



-2H\ 



(9) 



where the Hubble parameter H = H((f) is viewed as a 
function of <f>, V(4>) denotes the inflaton potential and a 
prime and dot denote differentiation with respect to the 
inflaton field and cosmic time, respectively. The energy 
density of the universe is p(cf>) = 3H 2 ((/>). It is assumed 
implicitly and without loss of generality that the inflaton 
varies monotonically with cosmic time such that <f> > 
(H' < 0). 

The Hubble slow-roll parameters are defined by 



H' 2 



''IT' 



and the spectral index is given by 



(10) 

where higher-order derivative terms in J\f = — \n(aH) 
have been neglected. To lowest-order in the slow-roll 
approximation, 1 — n s = 4e — 2n. This condition is a 
second-order, non-linear ODE for the dependent variable 
H = H(<P): 



(1-ef 



(I - e 2 ) fdlne 



[E1 
H 



12 



8 



H 2 



(1 



(12) 



Considerable insight into the nature of the general so- 
lution to Eq. (|12[) may be gained by expressing the en- 
ergy density of the universe in terms of the gradient of a 
'potential' function, W = W((f)), such that 



H 2 (q>)=4H 2 W'(q>), 



(13) 



where Hq is an arbitrary constant. Substituting this def- 
inition into Eq. (|12l) yields 



SW4>)] 



W"> 

w 



W" 

w 



A 2 
2 : 



(14) 



where A 2 = 1 — n s . 

The left-hand side of Eq. (|14l) is the Schwarzian deriva- 
tive of the function W((f>) and it is interesting that it 
arises in this cosmological context. Consequently, we may 



3 



now determine from Eqs. (p~3|) and ([14]) the general forms 
of the Hubble parameters for the full family of inflation- 
ary cosmologies that generate a constant spectral index. 

Since current observational bounds on the spectral in- 
dex inferred from the WMAP7+H0 data set are 0.939 < 
n s < 0.987 at the 2a confidence limit [8J, we focus on the 
red perturbation spectrum, A 2 > 0. A particular solution 
to Eq. (jT4"]) is W{4>) = exp(A</>) and the general solution 
is therefore given directly by 



w(4>) 



(IX: 



A0 + 5 



^ + d 



4Ai/ 2 



,A0 



(ce x <t> + dy 



(15) 



Moreover, comparing Eq. (|T4l with Eq. (|7J) and Eq. 
(TIB"]) with Eq. ([5]) immediately implies that H 2 (<fi) satis- 
fies the KdV equation 

- (1 - n s )H 2 ' + H 2 '" + ^H 2 H 2 ' = (16) 



if the general solution to Eq. (1141) is restricted to satisfy 
the Riccati equation 



W = —. — w 2 . 

4 



(17) 



It may be verified that condition (fl7|) is satisfied if 



_ 1 _2a 
cd c 



ad = -be = 1/2, 
As a result, the solution (|15l) can be expressed as 



H 2 W 
for cd > and 

H 2 W = 



Hi X 2 sech 2 



-Hi A 2 cosech 2 





(18) 



(19) 



(20) 



if cd < 0, where we have specified c = \d\ without loss 
of generality and have restored the dependence on the 
Planck mass for future reference. (If \d\, Eqs. (flU|) - 
(|20[) can be recovered by performing a linear translation 
A -1 In \d/c\ on the value of the inflaton field.) 

Given the general form of the Hubble parameter, the 
inflationary potentials can be deduced directly from the 
Hamilton- Jacobi equation ©. It is straightforward to 
verify that these potentials correspond precisely to those 
derived in to lowest-order in slow-roll. From the cos- 
mic dynamical systems point of view, the late-time at- 
tractor of these models is the power-law solution e = A 2 . 
There are two different potentials depending on whether 
the initial value of e is greater or less than A 2 @ . Within 
the context of the present discussion, the power-law 
model is the seed solution W(4>) = exp(A</>) and the two 
models are characterized by sgn(cd), i.e., by sgn(H // ). 



Eqs. (|T9[) and (|20|) are both wave solutions to the KdV 
equation (|16p and the former has precisely the form of 
the non-singular KdV soliton ([2]). This suggests a direct 
analogy can be drawn between such a wave and infla- 
tionary cosmology. In such an analogy, the inflaton field 
plays the role of a characteristic, wave-like coordinate on 
a two-dimensional spacetime {x, t}, the speed of the soli- 
ton is determined by the deviation of the spectral index 
away from the scale-invariant, Harrison-Zel'dovich spec- 
trum n s = 1, and the amplitude profile of the soliton is 
parametrized by the energy density of the universe (in 
appropriate units). 

The first slow-roll parameter for the 'soliton' solution 
(fl"9|). is given by 



\2 

= ytanh 2 (X(f>/2) 



(21) 



and determines the tensor-scalar ratio, r = V^/Vg = 
16e, where V\ = -ff 2 /(647r 4 e) and V\ = H 2 /(<iTr 4 ml) 
are the amplitudes of the density and gravitational wave 
perturbations @, 0]- This parameter is bounded from 
above such that r < 8(1 — n s ), which is consistent with 
the current observational 2a bound r < 0.24 M. 



B. The Cyclic Universe 

During inflation, quantum fluctuations in the inflaton 
field become frozen on super-Hubble radius scales be- 
cause the comoving Hubble scale decreases with time due 
to the rapid, accelerated expansion of the universe. How- 
ever, the comoving Hubble radius can also decrease if the 
universe undergoes a phase of slow, decelerated contrac- 
tion driven by a negative scalar field potential. This is 
the basis of the cyclic universe scenario. As shown in 
[Tfj| . the expression (|11[) for the spectral index is invari- 
ant under the duality e — >• 1/e and this implies there 
exists a one-to-one correspondence between inflationary 
and cyclic models that generate identical spectral indices. 

This suggests that a similar gravity /soliton analogy 
may be established for the cyclic universe scenario. In- 
deed, in the Hamilton- Jacobi formalism of the cosmologi- 
cal Friedmann equations (|9|) , the definition of the Hubble 
parameter H = a/a in terms of the scale factor a implies 
that 



a'H' = -~aH. 
2 



(22) 



Integrating Eq. (|22|) then yields the dependence of the 
scale factor on the scalar field in terms of the quadrature 



a{4>) — ao exp 



i r 4 




~2J 


V H> 



(23) 



where ag is an arbitrary constant. 

However, Eq. (|2"2"j) is invariant under the simultaneous 
interchange H(<fr) ■(-> a(4>). If we therefore consider a 



4 



'dual' cosmology where the Hubble parameter is given 
by H{4>) — a{(j>), Eq. (1231) implies that the dual scale 
factor is given by the quadrature 



a{4>) = do exp 



i r" 




~2J 


a' 



(24) 



Since the seed cosmology {H(<f), a((f>)} itself satisfies Eq. 
(p2| . the duality between the two scenarios is given by 
the simultaneous interchange of the Hubble parameters 
and scale factors of the two scenarios when all parameters 
are expressed as functions of the scalar field [Tl{ : 



a(4>) = H(4>), H{4>) = a(<t>). 



(25) 



It is straightforward to verify that under this duality, the 
Hubble slow-roll parameter indeed transforms as e = 1/e. 

As a result of this duality, the analysis of Section 3.1 
applies directly to the cyclic universe. We conclude, 
therefore, that there exists a one-to-one correspondence 
between solutions to the KdV equation and the respec- 
tive cyclic cosmological models when the spectral index 
is constant. For the cyclic model that is dual to Eq. (TT9"1) . 
the speed of the soliton is once more determined by the 
spectral index, whereas the amplitude of the soliton is 
now proportional to the square of the cosmic scale fac- 
tor. 



IV. MODIFIED COSMOLOGY 

In recent years, considerable interest has focused on 
cosmological dynamics arising from modified gravity the- 
ories. This is motivated in part by open questions in 
early universe cosmology, such as the singularity prob- 
lem, as well as providing alternative scenarios to dark 
energy models. At a phcnomenological level, such modi- 
fications can be quantified by altering the standard form 
of the Friedmann equation such that 



3H 2 = pL 2 (p), 



(26) 



where L = L(p) is a given function of the energy density 
and is determined by the specific model in question. 

Of particular interest in such scenarios are scaling (at- 
tractor) solutions, since these enable the generic asymp- 
totic behaviour of a cosmological model to be better 
understood. Scaling solutions are characterized by the 
property that the energy densities of the component mat- 
ter fields scale at the same rate as the universe expands 
(contracts). Such solutions were classified in Ref. [l2| 
under the assumption that the matter content of the uni- 
verse is comprised of a self-interacting scalar field with 
potential V(cf>) and a barotropic fluid with an adiabatic 
index 7. By introducing a parameter 



A 



L V 



(27) 



it was found that A — constant is a necessary condition 
for a scaling solution to exist. In that case, there exists 
an attractor solution for A 2 > 37, where the relative con- 
tribution of the scalar field energy density to the total 
density of the universe is = 37/ A 2 . There exists a 
second stable solution if A 2 < 6 where the scalar field 
dominates the fluid. It was further shown that these so- 
lutions exist if the condition 



p" 



1 



dHL(p)} 
dp 



(28) 



is satisfied |12j |. 

In this Section, we investigate the conditions that allow 
for scaling solutions in modified gravity to be interpreted 
as KdV-type solitons. We proceed by rewriting Eq. 
in the form 



P_ 

P 



3p^ 
2 p 2 



P 



12 r 



dlnL 



dp 



Defining a new variable a = a{<j)): 
p = Apod', 



(29) 



(30) 



where po is an arbitrary constant, then transforms Eq. 
into a Schwarzian differential equation: 



SM4)] = ( v 



dhiL 



P- 



dp 



(31) 



where the dependence of the square bracket on a is im- 
plicit. 

We now look for particular solutions to Eq. (|3"Tj) that 
satisfy the condition that the Schwarzian of a is constant: 



S[a(4>)] 



A 2 
2 ' 



(32) 



The discussion of Section 3 implies that the energy den- 
sity of the universe will satisfy the KdV equation if 
a' = (A 2 /4) — a 2 . As we saw above, the are two possible 
solutions when A 2 > 0, depending on whether a' > or 
a' < 0. (The product poa' is assumed implicitly to be 
positive-definite.) 

Eqs. d3TJ) and (J32J) imply that 



dlnL 



P- 



dp 



2 ■ 



Since a 1 ' 



-2aa', it follows that 



A 



£_ 

Po' 



(33) 



(34) 



Consequently, Eq. (|33|) simplifies to the first-order ODE 
dlnL 1 1 



dp 



2 Po X 2 -p 



(35) 



5 



and integrating yields the solution 



L= 1- 



PoA 2 



1/2 



(36) 



where we have chosen the appropriate branch of the gen- 
eral solution and the constant of integration to ensure 
that the standard relativistic cosmology is recovered in 
the low-energy limit p <C pqX 2 . The modified Friedmann 
equation (|26[) is therefore given by 



3H Z 



1 



P 



PoA 2 



(37) 



Finally, the results of Section 3 can be carried over to 
deduce that the corresponding scaling solutions, when 
expressed in terms of the energy density, are given by 
the KdV wave-solutions 



p = poA 2 sech 2 (A(/>/2) 
when a' > (po > 0) and 

p = -p A 2 cosech 2 (A^/2) 



(38) 



(39) 



when a' < (p Q < 0). 

The Friedmann equation (|3"7T) arises in a number of cos- 
mological models that are directly motivated by quantum 
gravity considerations. When po < 0, the model cor- 
responds to the Randall-Sundrum braneworld scenario 
[13j , where our observable universe is interpreted as a co- 
dimension one-brane embedded in five-dimensional, Z2 
symmetric anti-de Sitter space. The coefficient poA 2 is 
determined by the tension of the brane [l3j . On the 
other hand, the case where po > results in the Fried- 
mann equation for the Shtanov-Sahni (S-S) bouncing 



braneworld 14| . In this model, the universe is again in- 



terpreted as a one-brane embedded in a five-dimensional 
spacetime sourced by a bulk cosmological constant, but 
the extra fifth dimension is now assumed to be timelike. 

Furthermore, the Friedmann equation (|37[) arises 
generically in loop quantum cosmological (LQC) scenar- 
ios when po > [15|. In LQC models, the quadratic 
corrections in the energy density originate from non- 
perturbative quantum geometric corrections and become 
important at high energy scales. Indeed, in such a frame- 
work, poA 2 determines a critical density 



PO^ — Petit 



V3 



167T 2 7 3 



PP1, 



(40) 



where ppi is the Planck density and 7 « 0.2375 is the 
Immirzi parameter [l5l ]. 

In effect, therefore, by focusing on the KdV equation 
we have arrived at three different cosmological models 
that are all inspired by quantum gravity effects. The 
corresponding scaling solutions (|3"5)) and (|3"9"|) are those 
found previously in [H| . We may now interpret these 
solutions as wave solutions of the KdV equation. The 



non-singular soliton solution reflects the non-singular na- 
ture of the S-S and LQC scaling solutions. Due to the 
nature of the quantum corrections, the universe collapses 
from infinity (cj> — > —00 ), undergoes a non-singular 
bounce at p — p cr jt (0 = 0), and then expands to infinity 
{4> — > +00). From the point of view of a stationary labo- 
ratory 'observer', such dynamics would be equivalent to 
the time-dependence of the soliton amplitude as the wave 
propagates (modulo the appropriate relationship between 
cosmic and laboratory times). The cosmic bounce cor- 
responds to the passing of the peak of the wave. In the 
LQC scenario, the 'speed' of the wave is parametrised by 
the fractional energy density of the scalar field and the 
barotropic index of the fluid, A 2 — 3j/Q^,. 



V. DISCUSSION 

In the present work, it has been shown for the first 
time that the KdV equation arises in a number of im- 
portant cosmological scenarios, including the inflationary 
universe, the cyclic universe, loop quantum cosmology 
and braneworld cosmology. In each model, cosmological 
solutions can be reinterpreted as wave-like solutions to 
the KdV equation, and this allows for analogies to be 
drawn between cosmic dynamics and wave propagation. 

For example, in the inflationary scenario, we have 
found that the ODE determining the spectral index of the 
density perturbation spectrum generated during single- 
field, slow-roll inflation is closely connected to both the 
Schwarzian derivative operator and the KdV wave equa- 
tion. In principle, this allows for the full family of infla- 
tionary models that generate a constant spectral index to 
be classified in a very straightforward manner in terms 
of solutions to the KdV equation. 

In the region of observational parameter space r < 
8(1 — n s ), a formal analogy was established between the 
non-singular KdV soliton and the inflating universe. In 
such a correspondence, the scalar field plays the role of 
a wave-like coordinate, the speed of the soliton is deter- 
mined by the value of the spectral index and the ampli- 
tude of the soliton is parametrized by the energy density 
of the universe. Due to the inflationary /cyclic universe 
duality, similar conclusions hold for the simplest version 
of the cyclic universe scenario, although in this case the 
amplitude of the soliton is related to the scale factor of 
the universe. 

The general conditions for scaling solutions in a class 
of modified cosmological models sourced by a scalar field 
and a perfect fluid were considered. Requiring that the 
integral of the cosmic energy density (with respect to 
the scalar field) has a constant Schwarzian derivative led 
naturally to a modified Friedmann equation that arises 
generically in the Randall-Sundrum and Shtanov-Sahni 
braneworld models [l3l [lij and loop quantum cosmol- 
ogy scenarios [Hj]. In all cases, the corrections to the 
Friedmann equation are quadratic in the energy density. 
For such models, the scaling solutions may be interpreted 



6 



as wave solutions to the KdV equation, where the cos- 
mic energy density is again analogous to the soliton wave 
amplitude. 

Finally, it is worth remarking that the KdV equation 
is closely related to the non-linear Schrodinger (NLS) 
equation. This equation admits solutions that determine 
the propagation of solitons in Bose-Einstein condensates 
[l7l fl8j . A Bose-Einstein condensate is the ground state 
of a gas of N interacting bosons trapped by an exter- 
nal potential. In the limit where the interaction between 
the atoms is sufficiently weak, the mean-field approxi- 
mation may be employed. In this case, the macroscopic 
wavefunction for the condensate, ip, satisfies the Gross- 
Pitaevskii (GP) equation (T^ |: 



ihdtip 



2m 



(41) 



where V(x, t) is the trapping potential and m is the mass 
of the atoms forming the condensate. The scattering co- 
efficient is given by g — 4irh 2 Na/m, where a is the (s- 
wave) scattering length. 

By employing sufficiently anisotropic trapping poten- 
tials, it is possible to reduce the condensate to a quasi- 
one-dimensional configuration. Typically, the potential is 
given by V(x) = il 2 x 2 /2, where the trap strength ft <C 1. 
To a first-approximation, therefore, the potential can be 
ignored. In this limit, the condensate becomes homoge- 
neous and the GP equation is identical to the integrable 
non-linear Schrodinger equation: 



ihdt4> 



2m 



(42) 



The defocusing NLS equation (|42|) . where g > 0, sup- 
ports dark soliton solutions [17( ■ A dark soliton is an en- 
velope excitation characterized by a dip in the ambient 
density and a phase jump across the density minimum. 
Such solitons have been observed in a variety of Bose- 
Einstein condensates in recent years (for a recent review, 
see [2(|). The solution is given by [18[ 



h/f 



n 1 



l^dsl 2 ) 



ds 



sech 



V\ 2 {.r - rl i 



(43) 



where n is the background density, v is the speed of the 
soliton, (x — vt) is its position and c = ^Jngjm is the 



sound speed in the condensate. The spatial extent of 
the soliton is characterized by the healing length £ = 
h/ y/mng = l/y/Aima. In general, the speed of the soliton 
is bounded from above by the sound speed, v < c. 

The density profile of the dark soliton (|43p corresponds 
precisely to the energy density of the inflationary universe 
(fIU|) and the LQC scaling solution (|38p . On dimensional 
grounds, we can make the identification 



— < — > (x - vt)y/na (44) 
mp 

and view the scalar field as a wave-like coordinate. Mod- 
ulo a constant of proportionality, the spectral index may 
then be identified with the speed of the soliton: 



(45) 



whereas in the LQC model, the speed is proportional to 
the kinemetic parameters: 



J_ 



1 - 



(46) 



The deviation of the spectral index away from the scale- 
invariant perturbation spectrum is proportional to the 
speed of the soliton relative to the condensate sound 
speed. The maximal speed of the soliton is attained in 
the limit of the scale-invariant spectrum, n s — > 1. For the 
LQC scaling solution, the maximal speed corresponds to 
an equation of state p — —p for the fluid, which is the 
limit of a cosmological constant. 

The above analogies are not intended to be precise, 
but they do nonetheless suggest a new link between grav- 
itational and non-gravitational systems might be estab- 
lished through the KdV equation. It would be interesting 
to formalise such analogies further to establish kinematic 
correspondences between cosmology and condensed mat- 
ter physics. 
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